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Abstract 

Advances in synthetic methods have spawned an array of nanoparticles and bio-inspired 
molecules of diverse shapes and interaction geometries. Recent experiments indicate that such 
anisotropic particles exhibit a variety of 'nonclassical' self-assembly pathways, forming ordered 
assemblies via intermediates that do not share the architecture of the bulk material. Here we 
apply mean field theory to a prototypical model of interacting anisotropic particles, and find a 
clear thermodynamic impetus for nonclassical ordering in certain regimes of parameter space. In 
other parameter regimes, by contrast, assembly pathways are selected by dynamics. This ap- 
proach suggests a means of predicting when anisotropic particles might assemble in a manner more 
complicated than that assumed by classical nucleation theory. 
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I. INTRODUCTION 



Classical nucleation theory assumes the formation of ordered structures from similarly or- 
dered nuclei PQI2]- Mounting evidence, however, suggests that many molecular and nanoscale 
systems form ordered structures in more complicated ways, first associating as metastable, 
often amorphous aggregates. Such nonclassical crystallization pHSj has been observed in 
systems of spherical colloids [6]-|8] and the globular protein lysozyme [9h11J. as well as in 
numerous simulation studies [121 US] • 

Computational and theoretical work P^1 - [TB] reveals one set of circumstances in which 
particles bearing isotropic interactions assembly nonclassically: when attractions are made 
sufficiently short-ranged, the system's liquid-vapor critical point is submerged (in a density- 
temperature phase diagram) within the regime of solid-fluid coexistence. In what appears 
to be an immediate kinetic consequence of this thermodynamics, randomly dispersed com- 
ponents possessing short-ranged isotropic attractions, cooled below the liquid-vapor criti- 
cal temperature, tend to assemble into ordered solids only after forming transient liquid- 
like phases. However, most real components, from proteins to ions [17j to the plethora of 
recently-synthesized nanoparticles p3] , interact via anisotropic or 'patchy' attractions. Sim- 
ulation work [T9H23] reveals assembly pathways of such components to be in general richer 
than those of their isotropic counterparts. Further, experiments indicate that anisotropic 
proteins can crystallize via a metastable dense phase outside the liquid-vapor coexistence 
regime [21], an observation bolstered by recent simulations 

Two important ideas underpin our understanding of nonclassical assembly. The step 
rule of Ostwald [26] states that metastable precursors of the stable phase may appear if 
those precursors are closer in free energy to the parent phase than is the stable solid. The 
conjecture of Stranski and Totomanow (ST) [27J is the closely-related statement that the 
precursors that emerge are those confronted by the smallest free energy barriers to their 
nucleation. While these ideas receive broad support [13], recent evidence suggests that 
dynamical effects can invalidate the ST conjecture [28]. For all but one-component isotropic 
particles, then, it seems that there exists no simple physical picture that predicts when 
particles might assemble in a nonclassical fashion. 

Here we propose a step in this direction by considering a microscopic model prototypical 
of a collection of particles bearing isotropic and anisotropic interactions. In Sections [IT] 
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and in we introduce this model and use mean field theory to determine its phase behavior. 
We summarize this behavior in Fig. [TJ In Section IV we focus on thermodynamic states at 
which the solid phase is stable. We ask how the solid emerges if one begins with a well- 
mixed system and considers Langevin evolution in a free energy space of bulk 'density' and 
'structure' order parameters. We find that under some conditions there exists a free energetic 
driving force for assembly of the ordered solid phase via nonclassical pathways. In such cases 
the free energy surface local to the homogeneous fluid phase is stable in one 'direction' of 
order parameter space, and unstable in the other direction. Consequently, density and 
structure order parameters evolve sequentially, rather than simultaneously. There also exist 
thermodynamic states at which no such bias exists. In such cases, assembly pathways are 
determined principally by order parameter dynamics. We summarize these observations in 
Fig. [2] We conclude, in Section |V} by discussing an extension of this model in which the 
assembly of a solid phase is induced by the formation of a solid intermediate. This discussion 
is summarized in Fig. |3} 



II. MODEL 



We consider a collection of particles that live on the sites i G {1, A^} of a (i-dimensional 
hypercubic lattice. The presence or absence of a particle at site i is signaled by the occupancy 
variable taking the value 1 or 0, respectively. Particles bear unit orientation vectors Si, 
which, for simplicity, we assume to rotate in a plane [39]. We impose an energy function 
T-L = J2iLi (^i^ J2j ^ij ~ f^^i^ ' where j runs over the z = 2d nearest neighbors of i, and /i is a 
chemical potential. The dimensionality d profoundly affects the nature of fluctuations within 
the model, but at mean field level serves only to scale the strength of its pairwise interaction. 
We choose the pairwise interaction Uij to be a minimal representation of particles able to 
interact both isotropically and anisotropically: 

Uij = -riiUj (J + QSi ■ Sj) . (1) 

Here J is the strength of the isotropic interaction, and Q is the strength of the anisotropic 
interaction. This model is designed to describe vapor- and liquid-like phases of small and 
large occupancy number, respectively, in which particle orientations Si are disordered, and 
a ferromagnetic solid-like phase of large occupancy number in which particle orientations 
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show a high degree of order (a related coupled Ising-Heisenberg model possessing particle- 
vacancy symmetry was studied in [2S])- We next derive the free energy of this model in a 
mean field approximation. In such an approximation (see e.g. |3U]) the fluctuating variables 
at a given site feel only the thermal averages of variables at neighboring sites. The effective 
field at a given site is to this approximation T-Les = —n {Jp + QS ■ t + fi) = Ucs — pn. 
Here n and S are fluctuating variables, and we have introduced the collective density- and 
structure order parameters p = (n) and r = (nS), respectively. These order parameters 
serve to distinguish phases of low and high density, and phases in which particle orientations 
are disordered or mutually aligned. For future notational convenience we also introduce 
the Ising-like density variable = 2p — 1; we will use both (p and p. Thermal averages 
are defined self-consistently through the relation (A) = Tr(APcq), where the equilibrium 
measure Peq = q-^e-'^^'=« with q = Tre~^'^'=ff = 1 + 27re^("^^+'^)lo(/3(5|T|). Here I„ is the 
^th Qj-f^Qj- modified Bessel function of the first kind; (3 = 1/T (we adopt units such that 
fee = 1); and the trace Tr(-) = J2n=o i {^".i / ^'^ + ^n,o} (■) has been carried out by aligning 
r with Cj.. The effective Helmholtz free energy per site is then fes{p,^) = E — TS, where 
E = |(t/efT)-/ip and -TS = r(lnPeq) = -{'Hes)-T\nq. Thus fes{p,r) = -^{Ues)-T\iaq, 
or 

/efr(p,r) = ^-{Jp' + Qr') 

- Tln[l + e'^(-'^+'^)lo(/3Qr)] , (2) 

where r = |t| and p = p + Tln27r. We consider Eq. ^ to have been divided through 
by dimensions of temperature, and all parameters in that equation to have been de- 
dimensionalized accordingly. Equations of state for the density and structure order pa- 
rameters can be obtained by minimizing the free energy, and read 

and 

Ii(/3Qr) 



The expressions ([2])-(|4]) describe phases of vapor (low density, orientationally disordered: 
< 0,r = 0), liquid (high density, orientationally disordered: > 0, r = 0) and sohd (high 
density, orientationally ordered: > 0,r > 0). In the following section we derive the phase 
diagrams shown in Fig. [T] Readers not interested in the details of these calculations should 



focus on Section IV, in which we ask how the sohd phase emerges if it is stable and if we 
start from conditions of moderate density without orientational order. 



III. MODEL PHASE BEHAVIOR 

We first focus on the phase behavior of the model when either the isotropic interaction 
or the anisotropic interaction vanishes. For Q = we recover from g - ignoring field- 
independent terms and introducing K = J/A, ficoex = —^K and h = ^ (/U — /icocx) - the 
Ising model free energy /i(0) = — T \ncosh[(3 (Kcp + h)]. We recover from ^ the 
equation of state (p = tanh [P{K(f) + h)]. These expressions caricature the thermodynamics 
of the liquid- vapor phase transition [3l]. For K = 0, Eqs. ([2])-(|4]) describe, at /i = /icoex, a 
continuous phase transition in k = Q/A from a fluid phase having r = = to a solid phase 
whose order parameter scales near the critical point Kent = l^"^ as t^oi ~ (k — ^crit)"'^'^^- 

The phase diagram for general values of K and k (for T = 1) is shown in Fig. [Tj^a) 
(henceforth we focus on the case fi = /icoex)- It identifies a homogeneous fluid phase H 
(0 = = r); a regime of phase-separated (PS) liquid L (0 > 0, r = 0) and vapor V 
(0 < 0, r = 0); and a solid phase S (0 > 0, r > 0). The solid phase is described by Eq. ^ 
with p = Psoi(t) = Tlo(4/3Kr)/Ii(4/3/tr). The points {K,k) = (1,0) and (0, 1) are continuous 
critical points; Ci and C2 are lines of continuous critical points; and F (which abuts C2) is a 
line of first order phase transitions. The line M delimits the limit of fluid metastability. The 
equation of the union of the lines M and C2 is 2K = [f3 — 1/ k)^^ In (2/3/t — 1). It was found 
by equating derivatives with respect to r, at r = 0, of each side of Eq. Q (with p = Psoi(t)). 

Panels (b) and (c) of Fig. [T] show phase diagrams in the density-temperature plane for 
two choices of K and k. Panel (b) describes a case {K = 1.5, k, = 0.6) in which the 
solid phase becomes stable only well below the liquid- vapor critical point. Expansion about 
r = of ([2]) with p = Psoi{t) reveals the onset of r to be continuous with temperature 
(see inset), scaling below the solid phase critical temperature Tc = 1.08 (obtained from 
/3cK (1 + tanh [K {k,^^ — /3c)]) = 1) as Tsqi ~ (Tc — T)^^^. The density of the solid phase at 
the critical point is Psoi(t — ?■ 0) = (2k/3c)^^ ~ 0.90. A different scenario is seen in Fig[l|c): 
here the solid phase becomes viable above the liquid-vapor critical point (and stable with 
respect to the homogeneous fluid phase below T ^ 1.1) and the onset of r is now first order 
with K, (see inset). Cases (b) and (c) loosely resemble phase diagrams of Lennard- Jones 
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particles, with distinct vapor, liquid and solid phases; away from /i = /icocx (not shown) the 
phenomenology of this model is more akin to that of isotropic potentials of shorter range [T3] , 
where only one fluid phase is stable. 

IV. PATHWAYS OF ASSEMBLY OF THE SOLID PHASE 

With the phase behavior of the model established, we turn to the question of how the solid 
phase emerges if it is stable and if the system is prepared in the homogeneous fluid phase 
H (0 = = r). We imagine this latter phase, which is of moderate density and possesses 
no orientational order, to describe a well-mixed system. The thermodynamic driving force 
associated with evolution of the bulk phase from H to the solid is connected to the stability 
of the free energy surface, in the vicinity of H, in the (j)- and r-directions of order parameter 
space. These stabilities can be assessed by Taylor expansion of Eq. (|2]). Retaining only those 
terms required for thermodynamic stability (and ignoring field-independent terms) we find 

/eff(p,r) ^ ]^K{l-l3K)<f + 2K{l-l3K)T^ 

Recall that K = J/4, k = Q/4:, and </> = 2p — 1. The coefficients Cnm = 
{n\m\)~^ d^d!^fes{4','^)\(l>,T=o are positive constants (for K, k > 0). The signs of the co- 
efficients of the quadratic terms determine the stability of the fluid phase H. We see by 
inspection that the fluid is unstable to perturbations of density below a temperature Tp = K 
(recall that /3 = 1/T), and unstable to perturbations of structure r below a temperature 
Tt- = K. While Tp is the liquid-vapor critical temperature, T^. is not in general equal to the 
temperature at which the solid becomes stable. Ordering temperatures for specified model 
parameters are labeled in Fig. [l](b,c) [ID] . 

If the ordering temperatures Tp and are different, and if the assembly temperature T 
lies between them, then there exists a thermodynamic driving force along a preferred direc- 
tion of order parameter space, or, in other words, a thermodynamic impetus for nonclassical 
ordering. We can visualize the thermodynamically preferred assembly pathway by assuming 
evolution of the order parameters according to the equations 

P = -rp5p/efr(p,r) (6) 
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and 

f = -V,drU{p,r). (7) 

We assume the order parameter mobilities T p and Tj. to be constant, and we imagine them 
to be directly related to particles' translational- and rotational diffusion constants, respec- 
tively. It is likely that these approximations hold best in the case of one-component molecular 
crystallization. In general, order parameter mobilities will depend on the order parameters 
themselves, particularly whenever slow dynamics is encountered. Such is the case, for exam- 
ple, in models of systems undergoing gelation [32j or vitrification ^33j; in systems in which 
strong bonds are formed (e.g. in zeolite synthesis [M]); and in binary mixtures that exhibit 
slow inter-species mixing [281 EH] . The dynamics considered here neglects several other im- 
portant features of real systems, such as the effects of spatial diffusion, interfaces, and of 
assembly-impairing kinetic traps. Interfaces confer a surface tension between bulk phases, 
and can render order parameter mobilities anisotropic. In future work we will assess the 
extent to which the effects of surfaces on assembly can be captured by a Ginzburg-Landau 
expansion of the model defined by Eq. ([T]), and whether such expansions offer an alterna- 
tive microscopic route to 'phase field' models of crystallization (see e.g. [Ml EZ]). Here we 
focus on the simple dynamics of Eqs. (|6| and ([7|. We argue that this dynamics reveals, 
importantly, the thermodynamic preference for time-dependent evolution of bulk order. 

In Fig. [2]^a) we show Langevin pathways at T = 1 and T = 0.25 for model parameters of 
Fig. [l|b). Interpreted literally, the classical notion of assembly describes an approximately 
straight line trajectory between start- and end points in a phase space of (0, t). By con- 
trast, at the higher temperature the nonclassical 'density-structure' pathway is dominant, 
regardless of order parameter mobilities (pathways for Ft- = 1 and Ft- = 16 nearly super- 
pose), because the fluid H is stable to perturbations of structure but not of density. At 
the lower temperature the fluid is unstable in both directions of order parameter space, and 
both classical- and nonclassical pathways can be taken, depending upon order parameter 
mobilities. The density-structure pathway, characteristic of certain proteins' crystallization, 
owes its existence to the liquid-vapor critical point, as in the case of isotropic interactions. In 
panel (b) we show preferred pathways at T = 0.9 for the model parameters of Fig. [l](c). Here 
the nonclassical 'structure-density' pathway, characteristic of some melts [3], is preferred, 
though rapid evolution of p results in near-classical behavior. 
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V. INTERMEDIATE SOLID PHASES 



The density-structure pathway in our model is driven by the hquid- vapor critical point. 
However, recent work [211 [25] suggests that crystallization can be induced by assembly of 
a dense phase possessing some of the symmetries of the crystal even above the liquid-vapor 
critical temperature. To rationalize such behavior within the framework discussed here we 
can add to Eq. ([T| the nematic interaction term AUij = —Q2ninj cos (26* jj), where 9ij is the 
angle between neighboring particle orientations. The effective dimensionless Helmholtz free 
energy density for this augmented model is 

/eff (p, -r,^) = I [Jp" + Qr^ + Q2t^') 

- Tln(l + e^('^^+^)X(r,u;)), (8) 
where I{t,uj) = J^^'' ^^^^ e/3Q'^cose+/3Q2<^cos(26»)_ -^qj-q ^ ^ {n cos {26)) is a nematic order param- 
eter. From this free energy we find, via Taylor expansion, the ordering temperature for u to 
be = K2 = (52/4. The phase diagram for K = 0.5, k = 0.6, /t2 = 1 is shown in Fig.|3](a), la- 
beled with the ordering temperatures Tp, and T^; we focus on assembly at T = 0.9 (arrow). 
Here we observe a stable ferromagnetic solid phase Si (0i, Wi, ri) = (0. 99, 0.88, 0.95) having 
free energy density —1.1, and an unstable nematic solid phase S2 (02, 1^2, T2) = (0.91, 0.81, 0) 
of free energy density —0.6. In the absence of the nematic coupling k,2 the ferromagnetic 
solid (shown by line S'^ in (a)) is not viable at T = 0.9. When /t2 = 1 it becomes stable, but 
because T lies above and below we observe (Fig[3]^b,c)) assembly of the ferromagnetic 
phase Si via the unstable nematic phase S2, along the u — (f) — T pathway. Thus, assembly via 
a dense intermediate phase, whose symmetries are partially commensurate with the stable 
solid, occurs well above the liquid-vapor critical temperature. While different in detail, this 
behavior echoes the notion of assembly via metastable ordered intermediates considered in 
Ref. here it occurs because the free energy structure local to the homogeneous fluid 
phase favors assembly of the unstable solid phase S2, rather than its stable counterpart Si. 

VI. CONCLUSIONS 

We have used mean field theory to study two models prototypical of particles able to 
interact isotropically and anisotropically. While the approach considered here neglects im- 
portant effects of surfaces, molecular detail and thermal fluctuations, it reveals that complex 
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behavior can be driven by bulk free energy alone. We find that for a broad range of pa- 
rameters the free energy structures of these models favor assembly of stable solid phases via 
intermediate phases, either amorphous or ordered. For other parameter choices, by contrast, 
assembly pathways arc determined principally by dynamical considerations. One can ob- 
serve in such cases classical pathways along which intermediate phases resemble the stable 
phase. The work presented here suggests a simple microscopic framework within which to 
rationalize and predict the assembly pathways of anisotropic particles. 
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VIII. FIGURES 
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FIG. 1: Thermodynamic phase diagrams derived from Eq. ([2]). (a) In the space of varying 
isotropic- (K) and anisotropic (k) interaction strengths, we show regimes of stable homogeneous 
fluid H (moderate density, orientationally disordered); phase-separated (PS) liquid (L: high density, 
orientationally disordered) and vapor (V: low density, orientationally disordered); and solid S (high 
density, orientationally ordered). Cartoons depict the nature of these phases. Critical points and 



the nature of the lines F, M, Ci^2 are discussed in Section III (b,c) Phase diagrams in the density 
(p)-temperature (T) plane for model parameters such that the solid phase emerges below (b) and 
above (c) the liquid-vapor critical point. We expect nonclassical ordering (when the solid is stable) 
for temperatures between the ordering temperatures Tp and Tr (marked). The insets to (b) and 
(c) show the emergence of solid order r as a function of T to be continuous and discontinuous, 
respectively. Langevin trajectories at the three marked temperatures (red, green and blue arrows) 
are shown in Fig. [2j 



12 



(a) mT=l 



T = 0.25m (b) 



T = 0.91 






^ — rr=~4 






IT 




1 


i 







0.5 ■ 



1 



FIG. 2: Thermodynamically preferred assembly pathways derived from Langevin evolution on the 
free energy surface Eq. ([2]), with initial conditions (i;^(0), r(0)) = (10~^, 10^^). Order parameter 
mobilities Fp and are set to unity unless otherwise marked. Top: order parameters versus 
time; bottom: assembly pathways (blue) plotted atop free energy contours (gray) with time as a 
parameter. Panels (a) show trajectories at two temperatures under conditions used to generate 
Fig. [l|b). At the higher temperature (T = 1), the nonclassical 'density-structure' pathway is 
favored thermodynamically, because the fluid phase (</> = = r) is unstable to perturbations of 
density <f) but not to perturbations of structure r. Trajectories generated using structural mobilities 
Ft- = 1 (solid blue line) and F,- = 16 (dotted blue line) almost superpose. Cartoons depict the 
nature of three points along the trajectory. At the lower temperature (T = 0.25), by contrast, the 
fluid phase is unstable in both directions in order parameter space, and no thermodynamic bias 
for nonclassical ordering exists. The trajectory followed depends on order parameter mobilities, 
(b) Assembly at T = 0.9 under conditions used to generate Fig. [l]^c). Here the structure-density 
pathway is favored thermodynamically. 
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FIG. 3: Thermodynamics (a) and thermodynamically preferred assembly pathway (b,c) derived 
from Eq. Q, for K = 0.5, k = 0.6, K2 = 1. Panel (a) identifies two solid phases, a ferromagnetic 
phase Si and a nematic phase S2, in addition to the liquid- vapor coexistence curve. The assembly 
pathway shown in panels (b) (order parameters versus time) and (c) (parametric plot in order 
parameter space) is generated at T = 0.9 by Langevin evolution on the free energy hypersurface 
Eq. (|8|, starting from ((;^)(0), r(0), a;(0)) = 10^^(1,1,1), with equal order parameter mobilities. At 
this temperature the fluid phase H (0 = r = w = 0) is unstable to perturbations of nematic 
structure uj but not to perturbations of ferromagnetic structure r (because T lies below T^^ and 
above Tr] see panel (a)). Assembly of the ferromagnetic phase therefore occurs via the nematic 
phase. Cartoons depict the nature of the solid phases. 
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